Anisotropic pseudo-potential for polarized dilute quantum gases 
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Anisotropic pseudopotential relevant to collisions of two particles polarized by external field is 
rigorously derived and its properties are investigated. Such low-energy pseudopotential may be useful 
in describing collective properties of dilute quantum gases, such as molecules polarized by electric 
field or metastable 3 P2 atoms polarized by magnetic field. The pseudopotential is expressed in terms 
of reactance (K-) matrix and derivatives of Dirac delta-function. In most applications it may be 
represented as a sum of traditional spherically-symmetric contact term and anisotropic part. The 
former contribution may be parameterized by a generalized scattering length. The anisotropic part 
of pseudopotential may be characterized by off-diagonal scattering length for dipolar interactions 
and off-diagonal scattering volume for quadrupolar interactions. Two-body matrix element of the 
pseudopotential in a basis of plane waves is also derived. 

PACS numbers: 03.75.Fi, 34.20.Gj 
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I. INTRODUCTION AND PROBLEM SETUP 

The concept of pseudo-potential, i.e. full interparti- 
cle interaction being replaced by some less complicated 
"effective" potential, plays an important role in many 
subfields of physics [Q. In particular, properties of tra- 
ditional Bose-Einstein condensates (BEC) may be well 
understood just in terms of a delta-function potential, 
with its strength determined by a single parameter - scat- 
tering length [0, 0] which characterizes low-energy scat- 
tering between two particles. A rigorous derivation of 
pseudo-potential for spherically-symmetric interactions 
has been carried out by Huang and Yang B. Here 
I extend their derivation to anisotropic interactions. I 
also evaluate a matrix element of the derived anisotropic 
pseudo-potential in the basis of plane waves; this matrix 
element may be useful in studies of many-body properties 
of quantum degenerate gasses. 

Unusual collective properties of bosons and fermions 
with anisotropic interactions have generated a consider- 
able interest over the last few years, see e.g. Refs. ||, 
|, |, H g [Uj [n], [[I H 0. Below I enumerate sev- 
eral systems where the results of my analysis may be 
applicable. First, Yi and You []5] considered an ap- 
plication of strong electric field to an atomic conden- 
sate. The electric field induces atomic electric dipoles 
and thus anisotropic dipole-dipole interactions between 
the atoms. Another novel systems where the anisotropic 
interactions dominate at large separations are heteronu- 
clear molecules ||. Here an application of electric field is 
required to freeze the rotations of the molecules and to 
align the intrinsic molecular dipole moments with the 
field. Magnetic dipole-dipole interactions are present 
even for well-studied alkali-metal atoms. These inter- 
actions may be amplified for more complex atoms like 
europium and chromium |l5|, |l6|, [l7| with larger mag- 



netic momenta of the ground atomic state. The influence 
of such magnetic dipolar interactions on the condensate 
properties was discussed in Ref. || . New systems where 
the anisotropy of interactions may be also of interest are 
metastable 3 i-2 alkaline-earth atoms placed in external 
magnetic field. Here the long-range forces are dominated 
by interactions of atomic quadrupoles p8| . ft should be 
noted that the application of external magnetic or elec- 
tric field is important in all these examples - the field fixes 
quantization axis and a condensate may be described in 
terms of a single order parameter. 

In all the enumerated examples the collision process 
may be formalized using Fig. Here we show a pair of 
identical particles interacting in the presence of external 
uniform field. The z-axis is chosen along the direction 
of the field and angle 9 determines orientation of colli- 
sion (interatomic) axis f with respect to the field. At 
large separations r the particles are polarized along the 
direction of the field. In the most general case, as a 
result of a collision, a change in polarization may oc- 
cur (e.g. dipole moment of a molecule could end up 
pointing in the direction opposite to the field). We will 
disregard these non-adiabatic collisions. Then the in- 
teraction between the particles may be described by a 
unique potential V(r,9). Without loss of generality, this 
axially-symmetric potential may be expanded into Leg- 
endre polynomials Pl(cos9) 



V (r, 9) = V sph (r) + V anis (r, 9) , 

cosy 



(1) 

(2) 



L=2,4, 



Here V sp h if) and V m i s (r, 9) are spherically-symmetric 
(L = 0) and non-spherical contributions respectively. Al- 
though all even L contribute, at sufficiently large r the 
anisotropic contribution may be dominated by a single 
L term. In particular, we will focus on two practically 
interesting cases - dipolar (L = 2) 
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Knis (r, 9) -> Vdd 



C 3 



P 2 (cos 9) , r — > oo , (3) 



2 



and quadrupolar (L = 4) 
Knis {r,6)^V QQ 



P4 (cos ( 



(4) 



interactions. In the above equations constants Cl+i 
parameterize strengths of interactions and are propor- 
tional to the squares of respective multipole moments 
(e.g. molecular dipole moment). 
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FIG. 1: Geometry of collision process. At large separations 
colliding particles are polarized along external field. During 
the collision the particles are assumed to follow a unique adi- 
abatic potential. 



Realistic interaction potentials V(r) are singular at 
small interparticle separations and this singularity leads 
to well-known problems [fl9| in formulating perturbative 
expansion for many-body properties. In particular, ma- 
trix element of the interaction in basis of free-particle is 
divergent. To remedy this problem, the full interaction 
potential V is usually replaced by a pseudopotential V. 
A rigorous derivation of pseudo-potential for spherically- 
symmetric interactions has been carried out by Huang 
and Yang [|J . Here I extend their method to anisotropic 
interactions. 

Previously, for anisotropic dipolar interactions, Yi and 
You |j| proposed the following pseudopotential 



Vdd = gS(r) 



-| P 2 (cos ( 



(5) 



Here the first term is related to spherically-symmetric 
part of the full potential (||) and the second contribution 
is simply the long-range dipolar interaction (Q). This 
pseudopotential has been employed in a large number of 
studies of properties of BECs with dipolar interactions, 
see, e.g. g |, g §, | [Hj, |fij g|. Although straightfor- 
ward to work with in applications, the pseudo-potential 
(||) has certain shortcomings. For example, it is not valid 
in the vicinity of resonances. Pseudo-potential derived 
here resolves these shortcomings. 

The goal of this work is to consistently develop a 
pseudo-potential method for non-spherical interaction 
potentials. It will be required that two-body wavefunc- 
tions obtained with the pseudo-potential V and full orig- 
inal potential V to be equal at large interparticle separa- 
tions. The derivation of this pseudo-potential is carried 



out in Section |l| Certain properties of the derived pseu- 
dopotential are discussed in Section HI and we specialize 
the discussion to dipolar and quadrupolar interactions in 
Section |y|. Matrix element of the pseudo-potential in 
free-particle basis is evaluated in Section Finally, Ap- 
pendix contains a derivation of certain low-energy limits 
of K-matrix in the Born approximation. 



II. 



ANISOTROPIC PSEUDO-POTENTIAL 



We consider a solution of the Schrodinger equation for 
a relative motion of two particles interacting through a 
potential V(r) 



2fi 



(V 2 + fc 2 )$(r) = V(r)$(r) , 



(6) 



where /i is reduced mass of the pair and k is the rela- 
tive linear momentum. We assume that at sufficiently 
large separations r > r c , r 2 V (r) — * for any direction f . 
We also assume that the particles are contained in some 
large volume with characteristic size much larger than 
the extent of the potential r c . Some arbitrary boundary 
conditions may be imposed on the surface of the enclos- 
ing volume. At r > r c the wavefunction <E> (r) may be 
expanded in free-particle solutions 



$00 (r) 



I in 



{aim ji (kr) - I3im.ni (kr)) Yi rn (0, ip) , (7) 



where ji (kr) and ni (kr) are spherical Bessel and Neu- 
mann functions respectively and a\ m and /?/,„ are inte- 
gration constants. 

Following Ref. Q the pseudo-potential is determined 
by acting with |- (V 2 + k 2 ) on the asymptotic form (0) 



hi 

2/' 



(21- 



1)!!(7+1) 8(r) 



(8) 



Thus the original potential V(r) is replaced by a sum over 
"lumped" multipole sources placed at r = 0. To complete 
the construction of the pseudo-potential, we need to de- 
termine coefficients (3i m in terms of $ (r) . First we relate 
the integration constants ati m and (3i m by requiring the 
complete solution $ to be regular at r = 



Pirn — ^ K\m a l'm' 
I'm' 



(9) 



Here K\™ are the elements of the reactance (or K-) ma- 
trix used to parameterize multi-channel scattering [ p0[ . 
It is worth noting that the entire dependence of the pseu- 
dopotential on the original potential will be encapsulated 
in matrix elements of the K-matrix. 

Let us review some properties of the K-matrix. First 
it is related to more familiar scattering matrix S via S = 



3 



(1 + ifC)(l — iJC)^ 1 and further to transmission or T- 
matrix through S = 1 — iT. For low-energy collisions 
JC ~ — kT. K-matrix is real and symmetric 



(10) 



For identical bosons (fermions) only even (odd) partial 
waves need to be considered. The K-matrix is diagonal 
in I and m for spherically-symmetric potentials 



tan Si 



(11) 



pseudo-potential of Huang and Yang B is subsumed in 
this equation. Indeed, for spherically-symmetric interac- 
tions the K-matrix is diagonal in 1,1' and m,m' and is 
expressed in terms of conventional phase shifts, Eq. (|ll|) . 
Upon substitution of Eq. ( |Tl| ) into pseudopotential (|13|) 
we recover as a limiting case the result of Ref . . As the 
pseudo-potential of Huang and Yang Q the anisotropic 
pseudopotential is non-Hermitian and velocity depen- 
dent. 

The pseudo-potential ( |l3"| ) may be separated into 
spherically-symmetric V sp h and anisotropic V an is parts 



where Si is the phase shift for partial wave I. Compared 
to spherically-symmetric case, anisotropic potentials ad- 
ditionally couple different partial waves. For example, 
dipolar interactions, Eq. (|3|), couple s and d waves so 

that [/Cdd]qo 7^ an d quadrupolar interactions (||) cou- 
ple s and g waves. The scalar part V sph (r) of the poten- 
tial (||) assures non- vanishing /Cqq. In practice, K-matrix 
for non-spherical interactions may be found from a solu- 
tion of coupled radial equations pl| . In particular, it may 
be shown that for potentials parameterized by Eq. (|^), 
K\™ oc S mm i , i.e. the K-matrix is diagonal with respect 
to magnetic quantum numbers. In Appendix, I derive 
some elements of K-matrix in the Born approximation 
for dipolar and quadrupolar interaction. 

At this point we related the integration constants a; m 
and ftim via elements of K-matrix, Eq. (|^). Further, in 
the low-energy limit kr c <C 1, the integration constants 
ai m may be expressed in terms of (r) B 



Ctlr, 



1 1 

~¥T\~k l 



dr J 



21+1 



J YC rn m^oo{v) dn 



_ _ (12) 
Finally, combining equations (||)-([12|) we arrive at a gen- 
eralization of pseudo-potential V for anisotropic interac- 
tions 

v*(r)= -|E, rfrf l ra » (a&T'sw) , (13) 

where 



K 



d_ 

dr 

l m 



21' + 1 



/+ 1 / y; m ,(n)$(r,n) dn 



V = V E 



ph 



v a 



where 



and 



V sp h — — — 



M <H 

Im 



(k) v\Z 



Va n is , , 



M 



0>(J'T 



(16) 



(17) 



(18) 



In simpli fyin g the anisotropic part we used symmetry 
property ( |10[ ) of K-matrix. 

Let us focus first on the spherically-symmetric part of 
the pseudopotential and in the following section we will 
consider anisotropic part of the pseudopotential for dipo- 
lar and quadrupolar interactions. We expect that as in 
traditional BECs of dilute atomic gasses with spherically- 
symmetric interactions, the effect of V sp h on collective 
properties will be determined mainly by s-wave contri- 
.bution, i.e. / = 0, m = term in Eq. (J17|) 



J r=0 



V, 



P h 



h 2 AT 00 
n ^00 , 

M k 



,00 

^00 • 



It may be shown that for realistic potentials the following 
low-energy limit is finite 



a, s = -lim'" 00 



/C°° (k) 



(19) 



this quantity is a generalized scattering length. With this 
definition, the truncated pseudopotential reduces to 



(14,) 



V sph $ (r) « 4tt ^a ss S (r) |- (r <f (r)) , 



and M = 2/i is a mass of the collision partner. 



III. SOME PROPERTIES OF 
PSEUDOPOTENTIAL 

The derived anisotropic pseudopotential ( |l3| ) is one of 
the main results of this work. The spherically-symmetric 



r=( where we used § (r) = <5 (r) / (47rr 2 ) consistent with 
Ref. [0. Finally, for sufficiently slowly- varying wave- 
( lfii)iction, |dlog<i>/(ilogr| <C 1, we recover conventional 
contact pseudo-potential 



V, 



sph 



4-7T — a ss 8 (r) 
M 



(20) 



widely employed in studies of BECs. 

Having discussed spherically-symmetric part of the 
pseudo-potential, in the following section we consider 
anisotropic part of pseudopotential (|l^) for dipolar and 
quadrupolar interactions. 
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IV. DIPOLAR AND QUADRUPOLAR 
INTERACTIONS 



At this point we derived anisotropic pseudo-potential, 
Eq.(||). We separated the pseudopotential into 
spherically-symmetric and anisotropic contributions. We 
found that the spherically-symmetric contribution re- 
duces to familiar contact term (p0[), widely employed in 
studies of Bose condensates; the only modification being 
an introduction of generalized scattering length ([l9|) . In 
this section we focus on the anisotropic contribution and 
illustrate some of its properties for dipolar and quadrupo- 
lar interactions of identical bosons. These interactions 
were defined in the introductory Section as potentials 
that at large separations r are dominated by 



V^i M)^-fP 2 (cos< 



for dipolar interactions and 



c 

— Pi (cos( 



(21) 



(22) 



for quadrupolar interactions. 

Anisotropic interactions mix different partial waves 
(Im) and (I'm!) via off-d iagonal elements of K-matrix. 
From examining Eq. (A2) in Appendix, one may deter- 
mine that angular selection rules lead to coupling of s 
and d waves for dipolar interactions and s and g waves 
for quadrupolar interactions. In the following we assume 
that the dominant anisotropic effect on condensate prop- 
erties arises due to these particular couplings. Therefore, 



off-diagonal scattering length or volume by solving cor- 
responding scattering problem. 

Finally the total truncated pseudopotential is given by 



V 



A k2 

~M ass ^' 



V. 



DD.QQ 



(25) 



with 



K D n ? s $(r) 

Mr) 



\ dr ) 



■ J P 2 (cos 0) $ (r) dn 



9^P 2 (costf) [,|(r$(r))]_ } 
for dipolar interactions and 



(26) 



1 Ti 2 a 



2 7 M k 2 

,9 



[6 (r) [(£)" r 5 / P 4 (cos 6) $ (r) dSl\ + 
5(8!)^P 4 (cos#) [&(r#(r))] r J 



(27) 



for quadrupolar interactions. Quantities 5(r)/r n may be 
recognized as n derivatives of the Dirac delta-function. 

The constructed pseudo-potential depends on the rel- 
ative momentum k, i.e. the potential is velocity depen- 
dent. In practice, the velocity-dependence is most easily 
treated in momentum representation and in the next sec- 
tion we evaluate matrix element of the derived pseudo- 
potential in the basis of plane waves. 



* anis 



h z a: 20 

' M k 3 { 



V. 



QQ 



M 



fc 5 



- 2 °o°) 



For dipolar interactions it may be shown (Refs. [^2], |2^] 
and Appendix of this paper) that the following low- 
energy limit is finite 



,DD 



lim 



k 



(23) 



We will call this quantity off-diagonal scattering length. 
Similarly, for quadrupolar interactions we may introduce 
off-diagonal scattering volume 



ffl W = - lim 
9 fe^O 



k 3 



(24) 



In the Appendix we employ Born approximation and find 



C, 



DD „ 


1 M 


sd 








QQ r- 





It is worth noting that the above results are valid only 
away from resonances. In a general case one has to find 



V. MATRIX ELEMENT OF ANISOTROPIC 
PSEUDO-POTENTIAL IN FREE-PARTICLE 
BASIS 

While considering effects of two-particle interactions 
on properties of a quantum many-body system, one may 
require a matrix element of the derived pseudo-potential 
in free-particle (plane-wave) basis. We define this matrix 
element as 

V-(ki,k a> ki,l4)= (28) 
f dr dr'e- iki r e- lk2 r 'V (r - r') e *iV k = r ' . 

The pseudo-potential V, Eq. (|l3|), depends on the mo- 
mentum of relative motion of the interacting pair k = 
h (p — p'), where p and p' are momenta conjugated to r 
and r' respectively. To separate center of mass and rel- 
ative motions, we change the variables to R= (r + r') /2 
and r 12 = r — r'. With such a substitution 



1 



(27T) 



3 0ki+k' 2 ,ki+k 2 



M 



E ^n' (k') J dr 12 exp [-ik ■ r 12 ] v l f exp [ikf. r 12 ] 



IV mm 



5 



Here we introduced two relative momenta 



VI. CONCLUSION 



k = i (k x - k 2 ) and k'=~ (ki - 14) . (29) 

The delta function d-k'+k' ki+k 2 ensures conservation of 
the total linear momentum. Further we use partial-wave 
expansion 

exp [ik'- r 12 ] = 4^ £ i^j h (fc'ru) ij* mi (&') 1W (fia) 

iimi 

and arrive at 

y dr i2 exp [-ik ■ r 12 ] v(™' (r i2 ) exp [ik' • r i2 ] = 



(4tt) 2 i 



2 -J'-J 



2/ + 1 



(fc') k l Yjt m , (k')Yi 



Finally, the matrix element of the anisotropic pseudo- 
potential may be expressed in terms of relative momenta 
as 



v-(ki,k 2 ,ki,i4: 



dk'j+k^ki+ka X 

©^(kx-ka),^-]^) 



with 



?i 2 1 



s(k,k') = -'±-L y, i 1 ' 



■l'-l^lm 



IV mm' 

21 + 2 



-Y{! m , (k 1 ) Y lm (fc) . (30) 



v fc'/ 2Z 

Let us once again specialize the discussion to dipolar 



and quadrupolar interactions. As in Section IV we as- 
sume that the dominant anisotropic effect arises due to 
couplings of s and d partial waves for dipolar interactions 
and due to mixing of s and g waves for quadrupolar in- 
teractions. The corresponding truncated matrix clement 
(|30|) may be represented as 



" (k ' k ' )w i!f {a * 



^(k,k')i 



(31) 



with T replaced by 



T DD (k,k') 



,DD 



V5P 2 (cos0*/) + 



y/5 



I rigorously derived anisotropic pseudopotential aris- 
ing in the context of adiabatic collisions of two particles 
polarized by external field. Such low-energy pseudopo- 
tential may be useful in describing collective properties 
of dilute quantum gases, such as molecules polarized by 
an external electric field or metastable 3 P 2 atoms po- 
larized by magnetic field. The pseudopotential is given 
by Eq. (p~3|) . It is naturally expressed in terms of reac- 
tance (K-) matrix. The potential is non-Hermitian and 
velocity-dependent. It worth noting that in the deriva- 
tion I did not require the validity of the Born approxima- 
tion as in Ref. |||. Rather I followed method of Huang 
and Yang || and at large separations demanded the 
equality of solutions of two-body Schrodinger equation 
with a full original potential and with a pseudo-potential. 
Thus, compared to Eq. (|B|) by Yi and You ||, the derived 
two-body pseudopotential is expected to be also valid in 
a vicinity of scattering resonances. 

I argued that in most applications the pseudopotential 
may be represented as a sum of traditional spherically- 
symmetric contact term and anisotropic part, Eq. fl2q). 
The former contribution may be parameterized by a gen- 
eralized scattering length (|19|). We specialized discus- 
sion of the anisotropic part of pseudopotential to dipo- 
lar and quadrupolar interactions and found that it can 
be characterized by off-diagonal scattering length a^ d D , 
Eq. (|23|), for dipolar interactions and off-diagonal scatter- 
ing volume a^Q, Eq. (pij), for quadrupolar interactions. 
Although in a particular application these parameters 
should be determined from a solution of multi-channel 
scattering problem, I have derived a^ d D and o$~* in the 
Born approximations. Keeping in mind many-body ap- 
plications, I have also derived two-body matrix-element 
in the plane- wave basis, Eq. (|30|). Thus in this work I 
have rigorously derived anisotropic pseudo-potential for 
polarized dilute quantum gases and investigated its prop- 
erties. 
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.F QQ (k,k') = a?. Q I 3(k') 2 P±(cos9 k 



—k 2 
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P4 (cos Ok 



APPENDIX A: OFF-DIAGONAL SCATTERING 
LENGTH a° d D AND SCATTERING VOLUME a% Q 
IN THE BORN APPROXIMATION 



for quadrupolar interactions. In these expressions a ss is 
a generalized scattering length ( |l9| ) and a^J 3 and a^® 
are off-diagonal scattering length and volume defined by 
Eq. (p3|) and Eq. (p4l) respectively. 



lere we obtain expressions for reactance (K-) matrix 
in the Born approximation. Using the derived K-matrix 
we estimate off-diagonal scattering length a^ D and scat- 
tering volume for dipolar (DD) and quadrupolar 
(QQ) interactions. 
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The full solution of the Schrodinger equation 
be represented as 



$(r)=^^ 



! 111 



Ulm ( r ) 
kr 



may 



(Al) 



It can be shown that the radial functions ui m (r) satisfy 
the following system of coupled differential equations 



dr 2 



1(1 + 1) 



+ k } uim (r) = 



l'm!)u Vm > (r) 



(A2) 



with 
{lm\V\l'm!) (r) 



dft Y lm (Q) V (r) Y Vm , (n) . (A3) 



It is convenient to introduce regular and irregular solu- 
tions of homogeneous radial equations 



F x (kr) 
Gi (kr) 



krji (kr) , 
—kr ni (kr) 



and corresponding standing-wave Green's function 
9i (r, r') = - 



d (kr') Fi (kr) , r < r' 
d (kr) Fi (kr') , r > r' 



Using these definitio ns, s olutions to the system of inho- 
mogenous equations flA2| ) regular at r = may be repre- 
sented as 



uim (r) = ai m Fi (kr)- 



dr'gi (r, r') 



By comparing with Eq. (|?]) and Eq. we arrive at 
expression for elements of K-matrix in the Born approx- 
imation 



fC 



h Jo 

(A5) 

The Born approximation generally does not hold for 
low-energy atomic collisions, since realistic interactions 
are singular at small r. However, for dipolar interactions 
You and co-workers |2^, found numerically that away 
from resonances Born approximation works well for off- 
diagonal matrix elements. Keeping this observation in 
mind, below we derive off-diagonal scattering length and 
volume introduced in the main body of the paper. These 
parameters were defined as low-energy limits 



,DD 



lim 
lim 

fc^O 



[K D X(fc) 



k 3 



for dipolar (DD) and quadrupolar (QQ) interactions. In 
Section | the DD and QQ interactions were parameter- 
ized as 



Vdd = -f-P2(cos( 
C 

= —P4(C0S( 



(A4) 

where constants a; m are chosen to satisfy some bound- 
ary conditions. In the spirit of Born ap prox imation we 
may find solution of integral equations (A4) iteratively 
starting from 



2u^2 (lm\ v \ l ' m ') u iVrt> [ ^) th( r se definitions and Eq. JS5| ), we arrive at 



DD 

1 sd 



l sg 



1 M 

6V5 fi 2 



1 



1260 h 2 



c 3 



(A6) 
(A7) 



uim (r) w a hn Fi (kr) . 
In the lowest order in V at large r one obtains 



It is worth emphasizing that these results were derived in 
the Born approximation. In general case, to obtain pa- 
rameters entering anisotropic pseudo-potential (|l^) one 
r 2^ r°° nas to numerically -isolve the system of equations ([A2]), 

uim (r) — * aimFi (kr)+Gi (kr) ai m — — / dr'Fi (kr) ^Sp^df^iljn,l^l(&^ inity of resonances. 
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